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Abstract. Macroscopic Wigner islands present an interesting complementary approach to explore the prop- 
erties of two-dimensional confined particles systems. In this work, we characterize theoretically and exper- 
imentally the interaction between their basic components, viz., conducting spheres lying on the bottom 
electrode of a plane condenser. We show that the interaction energy can be approximately described by a 
decaying exponential as well as by a modified Bessel function of the second kind. In particular, this implies 
that the interactions in this system, whose characteristics are easily controllable, are the same as those 
between vortices in type-II superconductors. 

PACS. 41.20.Cv Electrostatics; Poisson and Laplace equations, boundary-value problems - 68.65.-k Low- 
dimensional, mesoscopic, and nanoscale systems: structure and nonelectronic properties 



1 Introduction 

The recent development of in-situ imaging techniques has 
allowed new experimental investigations of two-dimension- 
al mesoscopic devices consisting in small numbers of in- 
teracting confined particles. For instance, in type-II super- 
conductors, SQUID microscopy pQ or the multiple-small- 
tunnel junction method 2 have made possible the track- 
ing of the vortices. Even more recently, stable and metasta- 
ble vortex configurations in superconducting disks were 
observed using a Bitter decoration technique Many 
other systems are governed by the same physics: an inter- 
particle interaction, a confining potential, and possibly a 
thermal activation 

In parallel with the experiments, the simplicity of the 
input ingredients has allowed the development of numeri- 
cal simulations to determine both the equilibrium configu- 
rations and the dynamics of these systems. However, when 
the energy levels are very close, the uncertainties inherent 
in the numerical simulations make the confrontation with 
the experiments necessary. 

The difficulties in controlling exhaustive sets of ex- 
perimental parameters in real mesoscopic systems, such 
as the pinning or the interaction strengths, have led us 
to devise an analogue macroscopic experimental set-up, 
consisting in macroscopic Wigner islands, whose charac- 
teristics are easily tunable. As illustrated in Fig. ^ our 
macroscopic Wigner islands are constituted by millimct- 
ric stainless steel spheres (of radius R = 0.4 mm and 



mass m — 2.15 mg), sitting inside a horizontal plane con- 
denser of height h = 1.5 mm. The bottom electrode of the 
condenser is a doped silicon wafer, whereas the top one is 
a transparent conducting glass, thus allowing a direct cap- 
ture of the position of the spheres by means of a camera 
placed above the experimental device. A metallic frame 
of height h c = 1.5 mm, intercalated between the two elec- 
trodes and in electric contact with the bottom one, con- 
fines the spheres. Note that the confining frame could be 
set to a different potential than the bottom electrode, but, 
for simplicity, in the following we will focus on the equality 
case. When a potential difference Vq is applied to the con- 
denser, the spheres become charged, repel each other and 
spread throughout the whole available space, electrostat- 
ically confined by the outside frame. The voltage Vq goes 
typically from a few hundred Volts to a thousand, above 
which the condenser breaks down. By mechanically shak- 
ing the whole system, we can thermalize it: the spheres 
acquire a brownian motion [5j and obey Boltzmann statis- 
tics, in which the shaking amplitude plays the role of an 
effective temperature 

We then have a system in which the interaction be- 
tween the spheres, as well as the confining potential, can 
be easily adjusted. This device has been previously used 
to study small confined systems OH]; in such systems, 
as it has been shown in numerous numerical studies 9, 
I1UU11| . the equilibrium positions and the dynamics are 
highly dependent upon the interaction between the parti- 
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Fig. 1. Section of the experimentai device. Here the frame is 
circuiar, as it wili be used in section 0] 



can be made with sufficient accuracy, the interaction en- 
ergy is weli fitted by a decaying exponential as well as by 
a modified Bessel function Ko, both being controlled only 
by two parameters: their amplitude and their screening 
length. In section^ we determine a simple approximation 
of the confining potential in the case of a circular frame, by 
considering the frame as a hedge of spheres over which the 
intersphere interaction energy can be integrated. A correc- 
tion of the amplitude of the confining potential must be 
introduced to take into account the height difference be- 
tween the interacting spheres and the frame. This is done 
in sectional where the equilibrium configurations for up 
to 30 spheres and a circular frame are calculated. Com- 
parisons with the experimental data allow to adjust the 
amplitude of the confinement and to validate the model. 



cles. The same experimental set-up can be used to study 
pinning problems in two-dimensional elastic lattices. In- 
deed, defect-free lattices of a few thousand spheres can be 
obtained |12| . Again, the behavior of the system, in par- 
ticular far from the equilibrium, highly depends on the 
interaction potential. 

The comparison of our experimental equilibrium po- 
sitions with the calculated ones, for up to thirty parti- 
cles confined by a potential with axial symmetry, had 
led us to think that our interparticle interaction poten- 
tial was logarithmic, at least within the experimental dis- 
tance range \7\. The equilibrium configurations were also 
determined with an elliptic confinement [3]: the ground 
state configuration for 17 particles has been compared 
with the one for 17 vortices resulting from the minimiza- 
tion of the Ginzburg-Landau free energy in a mesoscopic 
type-II superconductor 1 of the same geometry [T5|. It has 
been suggested to the authors of the numerical study that 
the configuration that they presented was not the fun- 
damental state but the first excited level, as they finally 
confirmed |14|. Let us remind that the intervortices inter- 
action is described by the modified Bessel function of the 
second kind Ko |15|. which has the following asymptotic 
behaviors: 



*o(r/A) r ~„ -ln(r/A) 




r/A 



(1) 

(21 



All these indications on the nature of the interparticle 
interaction are indirect evidences and depend on already 
existing and studied interactions. Thus, they need to be 
confirmed and quantified. 

In this paper, we calculate numerically this interaction 
by considering the electrostatic problem of two spheres in 
electric contact with one of the electrode of an infinite 
plane condenser. Section [3 is dedicated to the numerical 
resolution of the problem. First, the Laplace equation is 
formally solved by means of a multipolar expansion. The 
numerically obtained electrostatic potential allows to cal- 
culate the electrostatic energy of the system, from which 
the interaction energy between two spheres is derived. In 
the range of distances in which the numerical calculations 



2 Interaction between two spheres 

2.1 Determination of the electrostatic potential 

We consider an infinite planar parallel plates condenser 
of thickness h; its lower plate is kept at zero potential 
and its upper plate is kept at the fixed potential Vo- On 
the lower plate, in electric contact with it, sit two con- 
ducting spheres of radius R < h/2; the centers of the 
two spheres are separated by the distance d > 2R. We 
introduce a Cartesian coordinate system r = (x, y, z) hav- 
ing the z-axis orthogonal to the plates of the condenser, 
with z — (resp. z — h) on the lower (resp. upper) 
plate; we choose the x and y-axis such that the centers 
of the two spheres are situated at x = ±d/2 and y = 0. 
Furthermore, we introduce two local Cartesian coordinate 
systems r_ = , z_) and r + = (x+, y+, z+), cen- 
tered, respectively, on the left and right sphere (see Fig-EJ , 
and their corresponding spherical systems of coordinates 
(r_,6L,0_) and (r+,6>+, </>+): 



x± = x =F d/2 = 



■■ r± sin 0± cos <fi± , (3) 
V± = y = r± sin 6± sin </>± , (4) 
z± — z — R = r± cos 9±. (5) 



To determine the interaction between the conducting 
spheres, we must solve the Laplace equation 



V A V = 0, 



(6) 



for the electric potential V(r), with the boundary condi- 
tions that V = on the lower plate and on the surface of 
the spheres, and V = Vq on the upper plate. To this aim, 
we begin to write the total electric potential V(r) as the 
sum of the potential Vq z/h of the empty condenser plus 
a perturbation due to the presence of the two spheres 



V(r) = Vo 



v(r) 



(7) 



The normalized perturbation v(r) is the solution of the 
Laplace equation that is regular in the space inside 



P. Galatola et al.: Determination of the interactions in confined macroscopic Wigner islands: theory and experiments 3 




Fig. 2. Geometry for the electrostatic interaction of two 
spheres inside a parallel plates condenser. 



the condenser and outside the spheres, and that satisfies 
the boundary conditions 



v(x, y, z = 0) = v(x, y, z = h) = 0, 



(8) 



«(r_ = = - z{ r-~W-\ ( 9 ) 

„ „ i \ z(r+ = R, 9+) 
v(r+ = R, 9 + ,<f>+) = - V ± - +; , (10) 



f(|r| — > oo) = 0, 



(11) 



where z(r_ = R, 0_) [resp. z(r + = R, 9 + )] is the z-coordi- 
nate of the point on the left (resp. right) sphere of lo- 
cal spherical coordinates (r_ = R, #_,</>_) [resp. (r + = 
R, 0+ , 4>+)] ; explicitly 



z(r± =R,6±) = R(l + cos 0±). 



(12) 



To determine the perturbation v(r), we start by con- 
sidering the solutions gim(r) of the Laplace equation ®, 
regular at infinity, obtained by separation of variables in 
the spherical coordinates (r_ , 0_ ,(/>_) : 



/^(cosg-)/^ 1 
fl<m(r) = — x cos(m^_), 



(13) 



where the P*? m (cos6L), with £ = 0, 1, 2, . . . and < m < £, 
are the associated Legendre functions of the first kind |16| : 



p < m (*) = L 2?/j-( 1 - ta 



V2 



(14) 



The elementary solutions l|13|l correspond to the usual 
terms of the expansion of the electrostatic potential in 
multipolar moments |17| : they are singular only at the 
center of the left sphere, r_ = 0. The factors h e+1 in 
their definition have been introduced to make them di- 
mensionless. Note also that we retained only the harmon- 
ics cos(m</>_), that satisfy the symmetry with respect to 



the (x, z)-plane of the present geometry. To construct, 
starting from the multipoles l|13|) , a suitable complete ba- 
sis of harmonic functions for our situation, we make a set 
of infinite mirror images of the potentials 113( 1 through the 
two plates of the condenser z = and z = h, such that the 
boundary conditions (JSJ are automatically satisfied, and 
we symmetrize the resulting potential through the (y, z)- 
plane, as required by our geometry. The resulting base 
functions are therefore 



oo 

9em(x, y,z + 2kh) 



-gtm(x, V, -z + 2kh) + gi m (-x, y,z + 2kh) 
-ge m (-x,y,-z + 2kh) ; (15) 

they are regular in the volume inside the condenser and 
outside the spheres, since their only singularities are at 
the centers of the spheres and at their successive mirror 
images through the two plates of the condenser, i.e., at 
(x = ±d/2, y = 0,z = ±R + 2kh), with k = 0, ±1, ±2, . . .. 
Using this set of base functions, we decompose the per- 
turbation v(r) as 



v(r) 



EE 

e=o m=0 



WmGlmfr). 



(16) 



The unknown coefficients vi m depend on the distance d 
between the two spheres and can be determined by impos- 
ing the remaining boundary conditions @ and <|10[) on the 
two spheres. Since the basis functions H15|) are symmetric 
with respect to the (y, z)-plane, only one of the two con- 
ditions has to be imposed. This is most efficiently done by 
projecting the boundary condition on the set of func- 
tions P(i m (cos#_) cos(m'0_), that represents a complete 
set of orthogonal functions — with the required symmetry 
with respect to the (x, z) plane — on the left sphere. The 
resulting boundary equations are therefore 



EE 

i=0 m=0 



r 




f sin 9 


-d0- / d(f>- 


Jo 


Jo 



G tm (r-=R,0-,<f>-) x 



Pe m (cos6L)cos(m>_) 
4nR 



h 
AttR 



if t' = m' = 
if £' = 1 and ml = 



(17) 



3// 

otherwise 



To numerically determine the coefficients ve m , we trun- 
cate the sum in the left-hand side of Eq. I|17|) to a finite 
I = £ m , and we evaluate Eq. (fT*7l) for m' — 0,1,..., £' 
and £' = 0,l,...,£ m . This gives a set of (£ m + l)(£ m + 2)/2 
linear equations in the (£ m + l)(£ m + 2)/2 unknowns vg rn . 
The integrals in Eq. {T7J|, that determine the coefficients of 
this linear system, are computed numerically by truncat- 
ing the mirror images sum in Eq. I(15J) to a finite |fc| = k m . 
By varying £ m and k m , we check for the convergence of 
the expansion. 
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Fig. 3. Contour plot of the electrostatic potential V(r) in the 
plane y = for two spheres of radius R — 0.2 h at the center to 
center distance d = 0.6 h. The spheres are indicated in white. 
Darker shadings correspond to lower values of V(r). 




Fig. 4. Same as figure|3]but in the plane z = R, i.e., the plane 
parallel to the plates of the condenser and passing through the 
centers of the spheres. 



2.2 Interaction energy 

At fixed potential difference Vq , the mechanical interaction 
energy F between the spheres is equal to the opposite of 
the electrostatic energy stored in the condenser 



F = ~QV , 



(18) 



where Q is the charge stored on the upper plate z = h. 
The latter, for our infinite geometry, is actually infinite. 
However, its variation AQ due to the introduction of the 
spheres is finite: this is the charge associated to the po- 
tential perturbation Vqv{y) in Eq. (0. To compute it, we 
make use of the reciprocity theorem, according to which, 
given two charge distributions p(r) and p'(r) that create 
the potentials V(r) and V'(r), respectively, 



p{r)V'{r)dr = / p'(r)V{r)dr 



(19) 



the volume integrals being performed over all the space |18| . 
As (p 1 , V) system, we consider a plane parallel condenser 
having its lower plate z = at zero potential and upper 
plate z = h at the potential Vq] its potential distribution 



is thus V = V z/h for < z < h, V = for z < 
and V' = Vq for z > h. As (p, V) system we consider the 
potential distribution having V = Vq v(r) for < z < h 
and V = otherwise; its associated charge distribution 
has three contributions: a surface charge Ao~ {x,y) on the 
lower plate z — 0, a surface charge Aa + (x,y) on the up- 
per plate z = h, and two point-like charge distributions 
centered on the two spheres, at (x = ±d/2, y = 0, z = R), 
associated to the corresponding multipolar singularities of 
the perturbation v(r). The surface integral of the charge 
distribution Aa+(x, y) is equal to the charge variation AQ 
due to the introduction of the spheres inside the con- 
denser: 

/ Aa+(x,y) dxdy = AQ. (20) 

Jz=h 

Now: 



p'(v)V(v) dr = Q, 



(21) 



since p' is non-zero only for z 
zero. On the other hand. 



and z — h, where V is 



p(r) V'(r) dr 



V Q [ Aa+(x,y)dxdy+^- [ p(r)zdr. 

Jz=h n JO R 

The last volume integral in this expression is performed in 
the region fin made of two arbitrary volumes encircling 
the two centers of the spheres, where p(r) is concentrated: 
it equals the z-component of the dipolar moment of the 
charge distribution p(r) around the origin r = 0. Using 
Poisson's equation p(r) = — eoV 2 [Vo«(r)], where eo is the 
vacuum dielectric constant, by successive integrations by 
parts and application of the Green theorem, we obtain 



/ p(v) zdr = e V f 



v(r)z ■ v — z 



dv 
dv 



dS, (23) 



where the surface integral is performed over two arbitrary 
surfaces dftn of outward normal v around the two centers 
of the spheres and z is the unit normal in the ^-direction. 
Taking for dflu. two spherical surfaces centered at (x = 
±d/2,y = 0, z = R) and using the decomposition Qlfi[l. 
one obtains 



/ p(r) zdr = 8-ireoVoh 2 



vqo- 



R 



Vio 



(24) 



In fact, the only terms in the multipoles (|13(l having a 
non-zero z-component of the dipole moment with respect 
to the point r = are the monopole (charge) term £ = 
m = and the dipole term £ = 1, m = 0. Finally, putting 
together Eqs. fH ty -ffiZ )! and the charge variation due 
to the introduction of the spheres at constant potential Vq 
can be expressed as 



AQ = -8ire V h 



R 



"00" 



"10 



(25) 



From Eq. (|18fl . the variation of the mechanical energy of 
the spheres is then 



AF(d)=e V 2 hf(d), 



(26) 
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Fig. 5. Normalized mechanical interaction energy Af as a 
function of the normalized distance d/h between the centers 
of the two spheres for 4 different radii of the spheres: R — 
0.1 h (dots), R — 0.2 h (squares), R = 0.3 h (triangles), R = 
0.4ft (diamonds). The lines indicate the exponential fit of the 
numerical data. 



with the normalized mechanical energy 



f(d) = 47T 



2.3 Numerical results 



R 

woo -r + "io 
h 



(27) 



In Figs. |3 and 0] we show typical cartographies of the elec- 
trostatic potential. As it is apparent from Fig. [3J above 
the top of each sphere the equipotentials are uniformly 
squeezed over the whole remaining thickness. Thus the 
wavelength of the corresponding perturbation along z is 
comparable with the thickness h of the condenser. Since 
V 2 V = 0, this implies a lateral relaxation on the same 
length scale: we thus expect that the interaction between 
the spheres is short ranged with a decay length « h (rather 
than R). 

In Fig. [3] we present the numerically computed me- 
chanical interaction energies Af(d) — f(d) — f(oo) as a 
function of the distance d between the two spheres, for 
different radii R of the spheres, from contact up to a dis- 
tance equal to 3.5 times the thickness of the condenser. 
Over this range of distances and for spheres of intermedi- 
ate size, roughly filling half of the height of the condenser 
(triangles in Fig. [SJ corresponding to R = 0.3 h), the nu- 
merical data are quite well fitted by an exponential law 



Af(d) = AeM-d/X), 



(28) 



with the decay length A ~ 0.29h. For smaller or larger 
spheres, this exponential behavior is violated at small dis- 
tances, while it is recovered at distances larger than h, 
with a decay length essentially independent of the size 
of the spheres. The amplitude A of the interaction grows 
with the radius of the spheres; as shown in Fig. it can 




Fig. 6. Amplitude A of the exponential fit of the mechanical 
interaction energy (see Ea. l28|l as a function of the normalized 
radius R/h of the two spheres. Dots: numerical data. Contin- 
uous line: power-like fit A = 852 (R/h) 5 - 5 . 



reasonably well approximated by the power-like behav- 
ior A ~ 852 (R/h) 5 - 5 . 

Nevertheless, we find that within the explored range 
for the distance d, which corresponds to the experimental 
situations, the numerical data are also well fitted by a Kq 
law 

Af(d) = A'K (d/\'), (29) 

as it is shown in Fig. [7| for the radius R — 0.267 h corre- 
sponding to the real radius of the spheres in the experi- 
mental device. The decay length is then slightly higher, 
A' ~ 0.32 h; the behavior of the corresponding ampli- 
tude A' as a function of the size of the spheres is shown 
in Fig. |H1 again, it is well fitted by a power-like behav- 
ior A' ~ 876 (R/h) 5 with essentially the same expo- 
nent. Indeed, the presence of a sphere inside the con- 
denser introduces a perturbation of the electrostatic po- 
tential that can be expanded in Fourier harmonics along z. 
The first harmonic varies along z as sinfaz/h): far from 
the axis of the sphere, the corresponding solution of the 
Laplace equation that has cylindrical symmetry along its 
axis and that decays to zero at infinity is proportional to 
sm(irz/h)Ko(Trd/h), where d is the distance from the axis 
of the particle. We thus expect A = h/ir ~ 0.318 h. 



3 Interaction with a confining ring 

In our experiment, the spheres are laterally confined by a 
frame that is in electric contact with the lower plate and 
almost touches the upper one. Here, we focus on the case 
of a circular ring of radius R c and height h c . To treat this 
confinement using the results of the preceding section, we 
model it as a necklace of touching spheres of diameter h c , 
whose centers lie on a circle of radius R c . We make the 
additional simplifying assumption that the electrostatic 
interactions are pairwise additive. These two hypotheses 
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Fig. 7. Normalized mechanical interaction energy Af as a 
function of the normalized distance d/h between the centers 
of the two spheres for the experimental radius R = 0.267 h of 
the spheres. Dots: numerical results. Continuous line: fit of the 
numerical data with the modified Bessel function Kq. 



A 1 = 




Fig. 8. Amplitude A' of the Ko fit of the mechanical inter- 
action energy (see Eq. 129ft as a function of the normalized ra- 
dius R/h of the two spheres. Dots: numerical data. Continuous 
line: power-like fit A' = 876 (R/h) 5A . 



ring. The screening length A, which does not depend on 
the radius of the spheres, is the same as the one deter- 
mined in subsection 12.31 The total normalized confining 
potential can then be written as 



Vo 



^exp 



n=l 




(31) 

where M = 2irR c / h c is the number of the necklace spheres 
making up the confining ring. To obtain a simpler analyt- 
ical expression, we approximate the sum by an integral. 
This is well justified in the limit R c 3> h c , which is already 
implicit in the hypothesis that the distance R c — r at the 
confining ring is large with respect to its height h c . In this 
limit, the confining potential becomes 



Mv Q 
2tt 



2tt 

exp 



-A" V 7 " 2 + R c - 2rR c cos< 



( 32 ) 

This integral cannot be expressed analytically in terms of 
elementary functions. However, a reasonably good approx- 
imation for r ~ R c /2 can be obtained by the saddle-point 
method JU] 



v c (r ~ R c /2) ~ MvoJ ^- ( - - ) exp 



2tt 



Rc 



Rr 



A 



_ (33) 

Finally, as R c 3> A and since it is found numerically in 
that case that the logarithm of the confining potential (|31|l 
has a dependence on r that is not far from linear, we 
approximate it with the first-order Taylor expansion of 
the logarithm of the approximate potential (|33|l around 
r = Rc/2. The resulting approximate confining potential 
that we shall use in the following is then 



v c (r) 



V2ir\R c 
vq exp 



2r 

f 

Rc 



r-Rc 



A 



(34) 



This expression fits very well the discrete sum <|3"T|) in the 
whole range of validity of the latter. Note that vo still 
needs to be determined, which will be done in the follow- 
ing section, where the approximations made above will be 
validated. 



are justified for distances to the confining ring large with 
respect to its height. 

For simplicity, as normalized interaction energy v c o be- 
tween a sphere and one of the fictitious spheres of the 
confining ring we take the exponential approximation 

v c o = w exp(-d/A), (30) 

where the amplitude vq depends on the radius of the sphere 
and on the height of the ring. The center to center dis- 
tance d can be expressed as d = (r 2 + R 2 — 2rR c cos 0) 1 / 2 , 
as a function of the distance r between the center of the 
sphere and the center of the confining ring, and the an- 
gle 4> between the two spheres, seen from the center of the 



4 Equilibrium configurations 

In order to validate our model, we consider the equilib- 
rium states of small Wigner islands. More precisely, we 
will successively focus on the general configurations, then 
on the precise positions of the spheres and finally on the 
energetic differences between the stable and metastable 
states. 

We use the approximate intcrparticlc and confining po- 
tentials l|28(l and l|34|l to determine numerically the equi- 
librium configurations of our macroscopic Wigner islands. 
The theoretical predictions are compared with the exper- 
imental observations, as it was done in Ref. [Jj for some 
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N 


Experiments 


Theory 


N 


Experiments 


Theory 








c 



1 £ 


1 fi 1 1 
l-O- 1 1 


i fi 1 1 

l-O- 1 1 





1-0 


1-0 


1 Q 


1 fi 1 9 
l-O- 1Z 


1 fi 1 9 


7 


1 fi 
1-0 


1-0 


on 


1 fi 1 3 
1-0-10 


1 fi 1 3 
l-O-lo 


Q 
O 


1 7 
1- I 


1 7 
1- f 


91 


1- (-lo 


17 13 
1- J -10 


9 


1-8 


1-8 


22 


1-7-14 


1-7-14 


10 


2-8 


2-8 


23 


2-8-13 


1-8-14 


11 


3-8 


3-8 


24 


2-8-14 


2-8-14 


12 


3-9 


3-9 


25 


3-8-14 


3-8-14 


13 


4-9 


4-9 


26 


3-9-14 


3-9-14 


14 


4-10 


4-10 


27 


3-9-15 


3-9-15 


15 


4-11 


4-11 


28 


3-9-16 


3-9-16 


16 


5-11 


5-11 


29 


4-9-16 


4-10-15 


17 


1-5-11 


1-5-11 


30 


4-9-17 


4-10-16 



Table 1. Experimental and theoretical ground state config- 
urations for 5 < N < 30 spheres of radius R = 0.4 mm in a 
condenser of height h = 1.5 mm confined inside a ring of radius 
R c = 5 mm. The cases were there is a discrepancy between the 
experimental observation and the theoretical prediction are set 
in boldface. 




model interaction potentials. We recall that our Wigner is- 
lands consist of 5 < N < 30 spheres of radius R = 0.4 mm 
contained in a condenser of height h = 1.5 mm held at 
the potential difference Vq — 700 V; the spheres are con- 
fined within a disc of radius R c = 5 mm and height h c 
practically coinciding with h. 

The theoretical stable and metastable equilibrium con- 
figurations are obtained by searching numerically for the 
local minima of the total, pairwise additive, normalized 
interaction energy 

H= J2 Af(|r 4 -!•,■!)+ J2 ^(M), (35) 

l<i<j<N l<i<N 

where is the 2D position of the i-th sphere with respect 
to the center of the ring. According to Eq. Il2(i[l . the in- 
teraction energy is normalized with respect to eoV 2 h. To 
minimize H , we employ a conjugate gradient method |20| . 
starting from many suitably chosen different initial con- 
ditions, in order to explore a significant portion of the 
complex energy landscape and find the various relative 
minima. 

The stable and metastable configurations form pat- 
terns roughly constituted by concentric shells on which 
the spheres are located. We shall refer to such patterns 
by means of the notation (A^-A^i-A^-. . .), where Ni is the 
number of spheres in the i-th shell from the center. 

The total interaction energy l|35|) depends on three pa- 
rameters: the decay length A, the amplitude of the inter- 
particle interaction A [see Eq. (|28|) ]. and the amplitude of 
the interaction with the ring vo [see Eq. i[34Jl]. We take the 
exponential approximation with A = 0.29ft. and A = 0.68 
as determined by a best fit of the numerical interaction 
energy between two equal spheres of radius R = 0.267/i. 
The remaining parameter vq, that depends on the height 
of the confining ring, is determined by adjusting the ra- 
dius of the ground state configuration for 5 spheres, that 
consists in a single shell of 5 particles, to the experimen- 
tal radius r — 2.25 mm. We thus obtain vq = 0.47. We 



Fig. 9. Ground state configuration for N — 20 spheres. The 
external circle is the confining ring, the small circles represent 
the experimental positions of the spheres, drawn to scale. The 
plus are the center of the theoretical positions of the spheres. 

note that, contrary to what one would expect on the ba- 
sis that h c > 2R (i.e., the necklace confining spheres are 
larger than the interacting spheres), we have vq < A. How- 
ever, it must be noted that, as we pointed out, the approx- 
imate confining interaction energies (|31(l and l|34|l are jus- 
tified only in the limit R c 3> h c , a condition that is not sat- 
isfied in our experimental situation, where R c /h c ~ 3. In 
that conditions, trying to determine more precisely by 
studying the interaction of two spheres of unequal heights 
would be pointless, since the approximation of the con- 
fining frame as a wedge of spheres would still remain. 
Nonetheless, as we shall see, our approximate interaction 
energies reasonably well account for almost all our exper- 
imental observations. 

4.1 Ground state configurations 

In Table ^ we report a comparison between the exper- 
imental and the theoretical ground state configurations 
for a number N of spheres ranging from 5 to 30. All but 
three of the ground state configurations are correctly pre- 
dicted. In Ref. 0, these same configurations were com- 
pared with predictions derived from other model interac- 
tion potentials: the best agreement was found for a loga- 
rithmic potential, that fails to predict the correct ground 
state only in 4 cases, the same 3 cases as our present model 
(N = 23, 29, 30), plus the N = 28 situation. 

As for the previously studied interactions, the observed 
discrepancies, that occur for dense packing of the spheres, 
are irrelevant, mainly because of the extreme smallness of 
the relative energy difference between the ground and the 
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Configuration 


-Rcxp [mm] 


Rth [mm] 


7?i og [mm] 


5 


2.25 




2.25 




2.25 


1-4 


2.44 




2.45 




2.52 


1-5 


2.48 




2.50 




2.76 


6 


2.35 




2.43 




2.52 



Table 2. Radii of the outer shell of 4 different stable (5 and 1- 
5) and metastable (1-4 and 6) configurations. R exp : experimen- 
tal values; 7i t h: theoretical values according to our model; iii og : 
theoretical value for a logarithmic potential with a parabolic 
confining. The two theoretical values were obtained by tuning 
the confining potential such as to reproduce the experimental 
observation corresponding to the first configuration. 



first excited state, as for the N = 23 and N = 29 cases. 
Moreover, for dense packing the hypothesis of pairwise 
interaction could be too strong and steric effects might 
make the experimental determination of the ground state 
tougher. Let's mention finally that another possible source 
of error might be the incorrect treatment of the confining 
potential for spheres too close to the ring. 



N | Excited state | E th [J] [ AE th [J] | AE cxp [J] 

18 TlTr! 4.1 x 1(T U 1.8 x KT 11 4.8 x lCT 11 

19 1-7-11 4.7 x 1(T 9 0.9 x 10 -10 1.3 x 10~ 10 

20 1-7-12 5.3 x 10" 9 0.3 x 10" 11 1.5 x 10" 11 



Table 3. Theoretical energies E t h of the first excited state 
for three different numbers N of spheres and theoretical AE t h 
and experimental AE axp differences of the energy of the excited 
state with respect to the ground state. 

Table reports a comparison between the experimen- 
tal and theoretical excitation energies for N = 18, 19, 
and 20 spheres. We also report the theoretical absolute 
values of the interaction energies of the first excited states. 
The agreement between the theoretical predictions and 
the experimental data is qualitatively good, considering 
that the relative differences AE/E 2 are extremely small. 
This last comparison finally validate our model. 

5 Conclusion 



4.2 Relative positions in the configurations 

In Fig. we compare the experimental positions of the 
spheres in the ground state configuration for N = 20 
spheres with the predicted one. We stress that these re- 
sults are not a fit of the experimental data, the only ad- 
justable parameter of our model, the relative amplitude 
of the confining potential vq appearing in Eq. (|34|l hav- 
ing been set once for all by comparison with the situa- 
tion for N = 5 spheres. Again, given that the condition 
of validity R c /h c 1 of our approximate confining po- 
tential is not satisfied, the agreement between the experi- 
mental data and our simple model is rather encouraging. 
The same good agreement is found for all the ground and 
excited state configurations. In particular, in Table [21 we 
compare the experimental and theoretical radii of the first 
two states for N = 5 and N = 6 spheres, along with the 
analytical prediction for a logarithmic potential with a 
parabolic confinement. This confirms that our theoretical 
interaction is closer to the real one than the previously 
predicted ones. 



Using a semi-analytical method, we have performed a di- 
rect calculation of the interaction between two conduct- 
ing spheres lying on the bottom electrode of a plane con- 
denser. We find that, within a significant range, the in- 
teraction energy can be described by a simple decaying 
exponential, as well as by a Ko function, both being gov- 
erned by a screening length that is typically a third of 
the condenser's height. On the basis of this interaction, 
our theoretical predictions for small Wigner islands con- 
stituted by up to a few tens of spheres are in good accor- 
dance with the experimental observations, thus validating 
our simple model. 

This study definitively completes the description of our 
system of interacting spheres, where the density, the inter- 
action amplitude and the temperature are now well-known 
and can easily be tuned. Thus, this system is a macro- 
scopic experimental model that easily allows to explore 
the properties of two-dimensional confined systems, such 
as vortices in mesoscopic type-II superconductors. In par- 
ticular, the understanding of the dynamics of the vortices 
could be enriched in a complementary way by correspond- 
ing studies on the macroscopic Wigner islands. 



4.3 Energy levels 

Finally, the energy differences between the ground and 
the first excited states have been explored. Experimen- 
tally, the possibility to submit the system, through a me- 
chanical shaking, to an effective temperature T that has 
already been calibrated allows us to explore these different 
states (jJJ. Neglecting the higher excited levels, that have 
a much larger energy, the difference AE between the en- 
ergy E2 of the first excited level and the energy E\ of the 
ground state is then obtained by measuring the ratio of the 
respective mean residence times < t 2 > and < t\ > |21j: 

<Ii> K e AE/ kB T_ (36) 
< T 2 > 
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